Measurement of Mutual Coulomb Dissociation in y/s^N = 130 GeV Au+Au collisions 
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We report on the first measurement of Mutual Coulomb Dissociation in heavy ion collisions. We 
employ forward calorimeters to measure neutron multiplicity at beam rapidity. The cross-section 
for simultaneous electromagnetic breakup of Au nuclei at a nucleon-nucleon center of mass energy 
y/SNN = 130 GeV is o~mcd = 3.67±0.26 barns which is comparable to the geometrical cross section. 
The ratio of the electromagnetic to the total cross section is in good agreement with calculations, 
as is the neutron multiplicity distribution. 

PACS numbers: 25.75.-q,29.27.-a,29.40.Vj,25.20.-x 



An interesting aspect of the Relativistic Heavy Ion Col- 
lider (RHIC) is the high rate of y-hadron collisions it 
produces. Photons from the highly Lorentz contracted 
electro-magnetic field produced by the heavy ions of one 
beam collide with the nuclei of the other beam (for an 
overview see jL|). For example, the flux of equivalent 
photons with energies 2 GeV < E 1 < 300 GeV in the 
target nucleus rest-frame corresponds to a y-nucleus lu- 
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20, 



X s 



With such a high flux, 



nuclear dissociation is highly probable. The calculated 
cross-section for Single Beam Dissociation of 95 barns 
limits the maximum beam lifetime at RHIC Q . 

We report the first measurement of Mutual Coulomb 
Dissociation (MCD) whereby both beams dissociate elec- 
tromagnetically. The calculation of this cross-section 
uses an extension of the Weizsacker- Williams formalism 
since it is dominated by 2nd order 2 photon exchange 
j| fl. The total cross section a to t, for the nuclear break 
up of both beams, due to either Coulomb or strong 
(hadronic) interactions, was calculated in H. There it 
was argued that uncertainties in the calculation of the 
two interactions partially cancel resulting in a theoret- 
ical error of 5%. One can exploit this to derive MCD 
and hadronic cross sections from the data presented here. 
The MCD process is also of interest because it can be 
used for luminosity monitoring at RHIC and at the Large 
Hadron Collider || |). Finally an understanding of the 
photon flux generated in peripheral heavy ion collisions 
is necessary to understand the emerging field of 77, 7A 
and 7 pomeron interactions (7| [|. 

The presented data are collected from the 3 experi- 
ments BRAHMS, PHENIX and PHOBOS at RHIC with 
Au beams at ^s^n — 130 GeV. MCD results in emis- 
sion of a few nucleons, dominantly neutrons, with small 
(few MeV) kinetic energy in the nucleus rest frame. In 
the laboratory frame they therefore have small angular 
spread with respect to the beam and close to 65 GeV 
energy. We measured neutron multiplicities close to the 
beam direction with Zero Degree Calorimeters (ZDCs) 



that are common to all RHIC experiments. The topologi- 
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are measured and compared to calculations |3j, In 
addition, Beam-Beam Counters which differ for each ex- 
periment were used to detect produced particles at large 
angles. 

The ZDCs are small transverse-area hadron calorime- 
ters with an angular acceptance of \9\ < 2 mrad with 
respect to the beam axis. They are located in the 
beam line behind the beam bending magnets at the four 
RHIC interaction regions and measure the energy of un- 
bound neutrons in small forward cones around the beam. 
Charged particles are deflected out of the ZDC accep- 
tance by the beam bending magnet leading to a mea- 
surement of neutral energy with very low background. 
The layout of the ZDCs is shown in Fig. |l] and further 
details can be found in ref. |J. The interaction trigger 
required an energy deposit above 12 GeV in the left and 
right ZDCs. The energy scale was determined from fits to 
the one and two neutron peaks, which are clearly seen in 
Fig. ^, but the detector response linearity was confirmed 
over the full range of observed energies. The energy res- 
olution at the one neutron peak is approximately 21%, 
consistent with test beam results ||. This agrees with 
our calculations which show that there is negligible en- 
ergy deposit other than from neutrons in the ZDC. 

The following discussion is based on the PHENIX anal- 
ysis, while the BRAHMS and PHOBOS results were ob- 
tained from similar analyses |10, 11 . Events were se- 



lected by a ZDC coincidence corresponding to a distance 
of < 20 cm from the nominal interaction point (calcu- 
lated using the time-of-flight to the ZDCs). This cut is 
identical to the selection used in ref. Jr| and 160,601 
events satisfied these requirements. 

The PHENIX Beam Beam Counters (BBCs)@ arc 
arrays of quartz Cerenkov detectors which measure rela- 
tivistic charged particles produced in cones around each 
beam (3.05 < \rj\ < 3.85, with 27r azimuthal cover- 
age). There are 64 photomultipliers (PMTs) in each BBC 
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FIG. 1: ZDC location: (a) top view (b) view along the beam 
axis. Nominal position and spread of neutrons, Au ions and 
protons at the ZDC location are shown. 

arm. A coincidence of the 2 BBC arms was used to tag 
hadronic collisions for the ZDC trigger sample. 

In Fig. ||(a) we plot the neutron multiplicity measured 
in the ZDCs (i.e. total amplitude relative to the single 
neutron peak) versus the BBC multiplicity (total ampli- 
tude relative to that of a single hit) Jl4j |. This amplitude 
was then rescaled to set the maximum multiplicity equal 
to 1 . We define 2 classes of events according to the num- 
ber of hit BBC photomultipliers above threshold, ubbc 

1. "Hadronic" events for which ubbc > 1 in each arm. 

2. "Coulomb" events with ubbc < 1 in at least one 
arm. 

The BBC trigger efficiencies for Coulomb and hadronic 
collisions were calculated from simulations of particle 
production for each class. We find that 99% of Coulomb 
interactions satisfy the Coulomb cut criteria since only 
3% of the events have a charged particle in either arm 
{§. However PHENIX calculated that only 92 ± 2% of 
hadronic interactions passed the "Hadronic" cut defined 
above [jl3| . The other events are misidentified and repre- 
sent a 14% contamination of the "Coulomb" sample. 

Diffractive hadronic interactions are not treated in our 
simulation and these are a potential background to our 
"Coulomb" sample. However this is a negligible correc- 
tion since, for example, the high energy p+W diffrac- 
tion dissociation cross-section is 20 mb |0|, which is a 
small fraction of the total p+W inelastic cross-section of 
2.7b. Background from accidental coincidences between 
single beam interactions was calculated from the individ- 
ual ZDC rates and found to be negligible. 

The ZDC multiplicity spectra are very different for the 
2 classes of events as can be seen from Fig. ||(b). The 
"Coulomb" events tend to have low total neutron multi- 
plicity whereas the fraction of "Hadronic" events with 1 
or 2 neutrons is very small. Hadronic interactions which 
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FIG. 2: (a) Scatter plot of total ZDC neutron multiplicity 
(sum of 2 arms) (measured as energy divided by the energy 
of one neutron) versus the BBC multiplicity, (b) Single arm 
ZDC neutron spectrum for the Coulomb and hadronic events. 



fail the "Hadronic" cut are peripheral events correspond- 
ing to a few elementary N-N collisions so events with low 
BBC multiplicity were used to model the shape of the 
14% background for "Coulomb" events. The corrected 
"Coulomb" shape is identical to the ZDC multiplicity 
spectrum obtained by measuring the left ZDC distribu- 
tion for events with one neutron in the right one. This 
cut strongly suppresses hadronic collisions. 

The ZDC trigger efficiency correction was obtained 
from the geometrical acceptance for neutrons in the ZDC. 
The angular distribution of neutrons about the beam 
direction is expected to be different for hadronic and 
Coulomb events so we discuss the two cases separately. 
In low energy photoproduction experiments the neutron 
multiplicity increased with ph oton energy and their max- 
imum pt was 120 MeV/c 0, [ItJ. For 65GeV neutrons 
this corresponds to a maximum opening angle of less 
than 2 mrad and therefore an efficiency of 100%. How- 
ever ~20% of the MCD cross section is due to equivalent 
photons with energy > 2 GeV for which no photopro- 
duction data are available. From the trend of the lower 
energy data we would expect that these events would 
have a mean neutron multiplicity larger than 1. Thus 
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FIG. 3: Left ZDC spectrum with "Coulomb" selection for 
events in which there is one neutron in the right ZDC. 



even if the px distribution broadens at higher photon en- 
ergies the ZDC inefficiency is expected to be small since 
it is unlikely that every neutron will miss the detector. 
Therefore, no ZDC acceptance correction is applied to 
the measured Coulomb rates. 

The ZDC acceptance correction for hadronic collisions 
was measured using an independent data sample which 
required only a BBC coincidence in the trigger. BBC 
trigger events were scanned for interactions in which the 
ZDC trigger condition failed. The measured efficiency is 
98 ± 2% and this correction is applied to the calculated 
hadronic rates. The acceptance correction was also in- 
dependently estimated assuming a neutron px spectrum 
characteristic of a Fermi distribution |l8j . The calculated 
geometrical acceptance per neutron in the ZDC detector 
is 75% at t i /snn=130 GeV. As can be seen from Fig. |^(a) 
multiple neutrons are detected for most nuclear collisions 
and so the calculated ZDC detection efficiency per event 
is close to unity, consistent with the measurement. 

We corrected our measured ratio of "Hadronic" events 
to total events to account for the (experiment dependent) 
BBC trigger efficiency, cbbc- For PHENIX ebbc wa s 
found to be [92 ± 2{syst)]% |l3j. The trigger efficiency 
is very weakly dependent on input assumptions of the 
Glauber model, such as the radius parameter. Note that 
it is only this efficiency and not the Glauber cross-section 
that is used in the current analysis. The BBC trigger had 
a background contamination of [1 ± l(syst)]% p3| . 

Table | presents our measured value of the ratio of 
Cgeom to the total ZDC cross section, a tot- This ra- 
tio agrees well with calculations in || [||. To compare 
to specific channels calculated for Mutual Coulomb Dis- 
sociation, we fit the measured ZDC energy spectra to 
neutron multiplicity distributions, convoluted with the 
energy resolution. We constrained the fits such that 
erE(2n) = y/2 x (7£;(ln), etc. This is justifiable from first 
principles. However the fitted areas can vary significantly 
if this constraint on the relative peak widths is removed 



a tot 



"tot 
"tot 

"(1,X) 
"(2,X) 
"(1,X) 



PHENIX 


PHOBOS 


BRAHMS 


a 


m 


- 


- 


- 


10.8 ±0.5 


11.2 


- 


- 


- 


7.1 


7.3 


.661 ±.014 


.658 ± .028 


.68 ± .06 


.67 


.659 


.117 ± .004 


.123 ± .011 


.121 ± .009 


.125 


.139 


.345 ± .012 


.341 ± .015 


.36 ± .02 


.329 




.345 ± .014 


.337 ± .015 


.35 ± .03 




.327 


.040 ± .002 


.042 ± .003 


.044 ± .004 


.041 ± .002 





TABLE I: Ratios of cross sections for experiment and theory. 
The values of a tot and u geom are in barns. 



and this variation dominates our systematic errors. 

Figure |^ shows the energy spectrum obtained in one 
ZDC (ZDC left) for the "Coulomb" event selection when 
the other ZDC pulse height is consistent with 1 neutron. 
The total number of events in Fig. 0, after background 
subtraction, corresponds to the cross-section for the (1,X) 
topology in which one neutron is observed in the right 
beam direction. Table | lists the sum of 1 neutron in left 
and right topologies as (1,X). Of course, any number of 
protons could also be emitted in the reaction but these 
are not detected. The decay topology (1,1) corresponds 
to exactly 1 neutron in ZDC left and 1 neutron in ZDC 
right. Similarly the fits to specific neutron multiplicities 
in the ZDCs such as (2,1) are derived from fits to the left 
ZDC spectrum, for events where we have cut on the pulse 
height in the right ZDC. The errors shown include both 
statistical and systematic errors in the fit procedure. 

For luminosity measurements the (1,1) topology is use- 
ful because it has a very low contamination from hadronic 
events, see Fig. |^(b), and only a small parameter depen- 
dence in the theoretical calculation |3|. Table | shows 
that for the ratio cr(l, l)/<Jtot the three experiments agree 
well with each other and the calculation. The error on 
the theoretical value is dominated by the error on a tot- 

In Fig. ^| we compare our measured 2n to In ratio to 
lower energy data on single beam dissociation fisj] . The 
ratio tends to increase with energy and the values are in 
reasonable agreement with calculations found in || 0]. 

Calculations of omcd suffer from the need to impose 
an arbitrary cutoff in the impact parameter bo below 
which electromagnetic processes are ignored and the in- 
teraction is assumed to be only hadronic. Reducing bo 
increases omcd but leaves a to t little changed. Also the 
error on atot is reduced due to the fact that at a given im- 
pact parameter the probability for either an electromag- 
netic or hadronic interaction cannot exceed one. Using 
this technique gives c ioi = 10.8 ± 0.5b at ^/snn = 130 
GeV [||. Since this result is more precise than current 
measurements of <r geom we use it to calculate (Jmcd, 

&mcd = &tot — cf geom = a to t ■ (1 — crg eom /cr tot ) (1) 



and the error on (Jmcd has contributions from the mea- 
sured ratios and the theoretical error on a to t- Using a 
weighted mean of the measurements of a geom /c7 tot by the 
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FIG. 4: Ratio of 2n to In cross-sections vs. ^/sjvjv. 

3 experiments in Eqn. [I] gives omcd = 3.67 ±0.26 barns. 

In Coulomb events each nucleus interacts indepen- 
dently with the field of the other. We therefore expect 
only a weak correlation between fragment multiplicities 
of the left and right going nuclei. For hadronic collisions 
both nuclei have the same average number of "wounded 
nucleons" and therefore the left and right energy in each 
event should be symmetric. We define the asymmetry 



A(ELeft, E Right) 



Pheft — E Right 
Eheft + E_ 



Right 



(2) 



where i?Le/t and ERight are the total ZDC energy signal 
in the Left and Right ZDCs respectively. Note we re- 
quire > 98GeV in at least 1 of the ZDCs to overcome the 
correlation inherent in the trigger. Figure || shows that 
hadronic events are much more symmetric than hadronic 
ones. We have also calculated the correlation function 
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PjELeft, ERjght) 

P{E Left ) ■ P(E R ight) 



(3) 



where PiE^^t, ERight) is the joint probability to have 
energies Eheft and ERight in each ZDC and P(ELeft) is 
the corresponding single probability. C2 is fiat and » 1 
for Coulomb events but reaches 4.6 for hadronic ones. 

In conclusion, we observe a class of events in Au+Au 
collisions at RHIC where neutrons are detected at beam 
rapidity but few particles are observed at large angles. 
The cross-section is comparable to the geometrical cross 
section and the ratio of the cross sections is in good agree- 
ment with earlier calculations of Mutual Coulomb Disso- 
ciation. Other features such as the neutron multiplicity 
distributions are also reasonably well reproduced. 
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